Measures of mixing quality in open flows with chaotic advection 
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We address the evaluation of mixing efficiency in experiments of chaotic mixing inside an open- 
flow channel. Since the open flow continuously brings new fluid into the limited mixing region, 
it is difficult to define relevant mixing indices, as fluid particles experience typically very different 
stretching and mixing histories. The repeated stretching and folding of a spot of dye leads to a 
persistent pattern. We propose that the normalized standard deviation of this characteristic pattern 
is a good measure of the mixing quahty of the flow. We discuss the link between this measure and 
mixing of continuously-injected dye, and investigate it using an idealized map. 



PACS numbers: 47.52. +j, 05.45.-a 

I. INTRODUCTION 

Mixing viscous fluids (whether Newtonian or not) in 
throughflow reactors is one of the most common and 
important industrial process [T]. Relevant quantitative 
measures of mixing efflciency are crucial for the design 
and evaluation of mixing systems. However, because of 
the wide spectrum of applications, there is not to date a 
universally accepted way of defining mixing efficiency [5] . 

Since the 1980s, mixing at low Reynolds number has 
been analyzed by fluid dynamicists in the framework of 
chaotic advection [21 [3], that is, the capacity of a mix- 
ing device to generate Lagrangian trajectories that sep- 
arate exponentially fast. In closed-flow mixing, that is 
when inhomogeneous fluid contained in a closed domain 
is stirred, some mixing indices stem directly from an un- 
derstanding of the mechanisms of chaotic advection [5]. 
Among classical measures, Poincare sections represent 
the region explored by a trajectory stroboscoped at dis- 
crete timesteps, and hence give the extent of the chaotic 
region [3H0] for time-periodic systems. Lyapunov expo- 
nents account for the stretching experienced by fluid par- 
ticles [IHS]. Topological entropy [T0Hl3] and braiding 
factors [T3] ITS^ describe the entanglement of Lagrangian 
trajectories. Such quantities, however, are hard to link 
directly to the degree of homogeneity achieved by the 
concentration field of a scalar to be mixed. More direct 
indices of homogeneity include the decaying variance of 
the scalar field [TBH^ or its entropy [H]. In some peri- 
odic flow configurations, the concentration field converges 
rapidly to a persistent pattern that repeats over time 
with decaying contrast P^^ES] . Such a pattern is a spe- 
cific eigenmode of the Floquet operator of the advection- 
diffusion equation, dubbed strange eigenmode; as such it 
allows an intrinsic characterization of the mixing flow. 

Less studied is the case of chaotic advection in open 
flows, in which fluid enters and then exits a limited mix- 
ing region. In this region, mechanical stirrers lead to 
exponential separation of fluid particles for the duration 
of their stay in the region (see Fig. [ija-b)). Extensive 
studies of idealized open flows in 2D have focused on the 



properties of the underlying dynamical systems, such as 
the multifractal set of never-escaping trajectories - called 
the chaotic saddle - and its influence on Lagrangian tra- 
jectories flowing through the mixing region |27H33j . Such 
studies focused mostly on reaction rates, and less on the 
mixing quality. 

On the other hand, a variety of practical indices of mix- 
ing quality for open reactors can be found in the chemical 
engineering literature [MH55] . These indices are mostly 
used for turbulent flows. The intensity of segregation was 
deflned by Danckwerts ^Sdj as the variance of the scalar 
field in the outflow. This index characterizes directly the 
homogeneity, hence the quality, of the end product. Re- 
fined measures of mixing take into account the spatial 
scale of the scalar field as well as its variance [HS]. One 
can also compute the departure of the residence-time dis- 
tribution (RTD) of fluid particles inside the mixer from 
a purely exponential RTD that characterizes a perfect 
mixer, in the language of chemical engineering. 

In recent work [39], we have studied the mixing dy- 
namics of quasi-2D free-surface channel flows stirred by 
two rods moving in an eggbeater protocol (see Fig. [ija) 
for a picture of the rods and the channel). A spot of 
dye released upstream of the stirrers (Fig. [ija)) is ad- 
verted by the main channel fiow and enters the mixing 
region where it is stretched and folded by the stirrers. 
Some dye particles leave the mixing region very quickly 
(Fig. [ijb)), while other parts of the spot remain in the 
mixing region longer (Fig.[ljc)). We demonstrated that, 
after a short time, the concentration field of the dye 
takes on a persistent pattern whose mean intensity de- 
cays with time '39' . By analogy with persistent patterns 
in closed fiows [22, |231 US] , we called such patterns open- 
flow strange eigenmodes. 

In the present work, we examine which information 
about mixing quality can be gleaned from such exper- 
iments in open fiows. In particular, we investigate 
whether relevant measures of mixing quality can be de- 
rived from the concentration field of the strange eigen- 
mode. This pattern is very robust, in the sense that 
it does not depend on the upstream location where dye 
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FIG. 1. (a) Experimental setup: a viscous fluid in a channel is 
stirred by two rods moving in an eggbeater motion. A spot of 
black dye is injected upstream of the rods, (b) Mixing pattern 
observed shortly after the dye has entered the vicinity of the 
rods: the rods stretch some dye particles into thin filaments, 
while others escape downstream without being stretched and 
mixed, (c) Mixing pattern observed at longer times: the in- 
tensity of dye in the pattern decays with time due to diffusion 
- this is the 'strange eigenmode' regime, (d) If the rods are 
rotated in the opposite direction to (a-c), they prevent fluid 
from flowing along the channel walls. Two stagnation points 
(plotted as disks) on the channel sides trap fluid in their vicin- 
ity for long times. In that case, the mixing pattern occupies 
the whole width of the channel. 



was injected. It is therefore natural to use the pattern to 
characterize the mixing device. 

The paper is organized as foUows. We first describe 
the experimental apparatus and methods in Sectionjll] In 
Section |III[ we describe the history of a spot of dye pass- 
ing through the mixer, and the evolution of the resulting 
concentration field. This helps to extract relevant infor- 
mation from the concentration images and time series, as 
we describe in Section |IV[ In particular, we introduce a 
new mixing index called the eigenmode index, which is 
the normalized standard deviation of the persistent con- 
centration pattern. In Section [Vj we relate this index 
to the intensity of segregation for a continuous injection 
of heterogeneity as introduced by Danckwerts [M] . The 
paper ends with a summary and discussion of the results. 



II. MIXING DEVICE CONFIGURATIONS 

Experimental data are obtained in an open-flow chan- 
nel whose experimental set-up has been described in more 



details in [321 EHj- Viscous fluid (cane sugar syrup) flows 
at a fixed rate through a long and shallow transparent 
channel. Two rod-stirrers (see Fig. [l]) travel on inter- 
secting circular paths at mid length of the channel, in an 
eggbeater protocol. Because of the low Reynolds num- 
ber (see below), the mixing flow is time-periodic. We 
select the eggbeater protocol for the simplicity of its ge- 
ometry as well as its ability to promote efficient chaotic 
advection. Indeed, the intersection of the trajectories 
of the rods ensures chaotic advection, due to stretch- 
ing and folding of fluids particles. Through the stir- 
ring frequency, we control the order of magnitude of the 
mean number of periods spent by fluid particles inside 
the mixing region, which was varied between 4 and 24 
in our experiments [391 . Depending on the stirring fre- 
quency, Reynolds numbers range from 2 to 10, well within 
the laminar regime. The mixing region may be defined 
loosely by the region spanned by the trajectories of the 
stirrers; a more rigorous definition will be given in the 
next section. 

The orientation of the trajectories of the rods (with 
respect to the direction of the main flow) determines the 
routes along which fluid flows preferentially inside the 
mixing region. We may rotate the rods so that they assist 
fluid in passing along the channel walls (Fig. [T|a)-(c)), 
where they drag fluid inside the mixing region only while 
on the upper and central part of their trajectory. Or 
we can rotate the rods in the opposite direction, so that 
fluid moving along the side walls cannot cross the mixing 
region, but is forced into a central funnel (Fig. flld)). 
We call the two cases respectively protocol type A and 
protocol type B. 

We perform mixing experiments by releasing a spot 
of dye upstream of the rods (Fig. [ija)). Pictures of the 
evolving dye pattern are taken at every period of the stir- 
ring motion with a high-resolution digital camera. Our 
set-up has been designed to measure quantitatively the 
concentration field from image processing (see [321 HO] 
for more details). This requires in particular using only 
backlighting of the dye pattern. 



III. HISTORY OF A DYE SPOT 

We now turn to a qualitative description of the differ- 
ent spacetime trajectories of fiuid particles, and of the 
resulting dye concentration fields. The understanding 
we gain here will be used in the next section for defining 
mixing indices. 

From the moment when the spot of dye is injected 
upstream of the rods, the fluid particles in the spot can 
have two types of radically different future histories: they 
either enter the stirring region and are mixed to some 
degree, or they pass through the region without being 
mixed. The latter behavior occurs predominantly in type 
A protocols, where fluid particles can flow freely along the 
channel walls. 

The patch of fluid particles that enter the mixing re- 
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gion at a given period consists of dyed particles as well 
as some undyed white particles. After a few stirring pe- 
riods, only white undyed fluid enters the mixing region. 
Before diffusion becomes important, the patch contain- 
ing dyed particles is stretched and folded by the rods. At 
each stirring period, the patch moves to a new part of 
the mixing region, and its previous location is replaced 
by fresh, undyed fluid that entered the region later. 

As time goes on, the original dyed patch inside the mix- 
ing region develops thin fllaments as it is stretched and 
folded by the rods (see the filaments inside the mixing re- 
gion in Fig. [ijb)). Nevertheless, as long as the width of 
the filaments is large enough, fiuid particles of the dyed 
patch are not mixed by diffusion with undyed fluid par- 
ticles. For example, in Fig. [ijb) most dye filaments bear 
the same concentration level as the initial dye spot. At 
every stirring period, some particles of the dyed patch 
leave the mixing region to compensate for fluid entering 
the mixing region; for short residence time they have not 
been mixed at all with undyed particles (see the dark fil- 
aments in the upper part of Fig.[l](b), downstream of the 
rods). 

For longer residence times, mixing by diffusion of fluid 
particles with different entry times occurs at locations 
where the patch has been compressed down to the diffu- 
sion scale at which stretching and diffusion balance - the 
Batchelor scale ^SJ HU |32] • The Batchelor scale is given 
by 

WB = V^M, (1) 

where k is the diffusion coefficient and A is the mean 
stretching (Lyapunov exponent). The scale wb is the 
smallest that can be observed in the mixing pattern, 
and increases with k. Once the patch has reached the 
scale wb , diffusion efficiently blurs its particles with par- 
ticles that have different entry times. In particular, 
dye filaments are blurred with surrounding white fluid, 
and intermediate gray concentration levels appear as in 
Fig. [l](c). Mixing is efficient for particles of the patch 
that have stayed inside the mixing region at least until 
their filaments reach w^. 

Such filaments approach the (non-elliptic) periodic 
points of the flow, which return to the same position 
after a given number of stirring periods and hence re- 
main forever inside the mixing region. The set of all 
periodic points is known as the chaotic saddle for open 
flows |5H1 Ell E2- It is a set of zero measure. The peri- 
odic points are typically located in between the regions 
occupied by the patch for short residence times, which 
cannot contain periodic points as they are filled with new 
fluid at each period. The mixing region can now be de- 
fined more precisely as the basin where the orbits of the 
chaotic saddle remain trapped forever, i.e., the convex 
hull of the set containing all these periodic orbits. Us- 
ing arguments from lobe dynamics [43], it can be shown 
that this set is delineated on the upstream side by the 
upstream boundary of the strange eigenmode pattern. 
On the downstream side, the extent of the mixing re- 



gion is bounded by the symmetric line of the upstream 
boundary with respect to the horizontal axis linking the 
rotation center of the rods. 

Dye filaments approach a periodic orbit of the chaotic 
saddle along the orbit's stable manifold - the set of fluid 
particles that converge to the periodic orbit. Conversely, 
they escape an orbit along its unstable manifold - the 
set of particles that diverge from the orbit. The typical 
situation is that fluid particles in the patch approach an 
orbit of the chaotic saddle along its stable manifold, then 
recede along its unstable manifold towards another orbit, 
and so on until they finally escape downstream. 

The dye pattern of Fig. [ijc) traces the unstable man- 
ifold of the chaotic saddle. As the unstable manifold is 
invariant under a full stirring period, dye always returns 
to the same pattern. We observe indeed in our experi- 
ments that patterns such as the one in Fig. [ijc) are per- 
sistent in time. Because of diffusion and the loss of dye 
downstream, the intensity of the pattern becomes weaker 
and weaker. 

In previous work |39j , we have shown that for protocol 
type A the dye pattern converges to a persistent con- 
centration field that repeats itself periodically and whose 
mean is decaying exponentially with time. The concen- 
tration field obeys the equation 

C(x, t^nT)^ exp{-nT/T) C(x) 

where t is the decay time (l/r is the decay rate) . (In gen- 
eral 1/t is complex, which leads to oscillations around a 
decaying trend; here we take its real part.) By analogy 
with closed flows [52], we call C'(x) the strange eigen- 
mode of the flow, as it is an eigenmode of the advection- 
diffusion operator [39]. Fig. [Ijc) displays an example of 
strange eigenmode. A physical picture of the support of 
the strange eigenmode (the space where |C'(x)| > 0) is 
given by the unstable manifold of the chaotic saddle - 
a fractal set of zero measure [31] - fattened by diffusion 
into filaments of width wb- Conversely, the complement 
of the support of the eigenmode does not contain any pe- 
riodic orbits and corresponds to the locations occupied by 
a patch of fluid shortly after its entry time, before it has 
reached the Batchelor scale. It can be seen in Fig.[2]that 
the short-time iterates of the dye spot are located inside 
the holes of the eigenmode. For example, the filamen- 
tary pattern in the mixing region in Fig. [2](a) is located 
in the holes of the long-time pattern of Fig. |2jc) . This 
is also true for the fluid escaping the mixing region: the 
escaping filaments in the lobes of Figs. [2]ja-b) are found 
within the holes of the pattern that periodically leaves 
the mixing region at long times (Fig. [2](c)). The prop- 
erties of the strange eigenmode were investigated more 
thoroughly in |39j . 

For protocol type B, however, we found that the pat- 
tern continues to evolve slightly with time and its decay 
does not exhibit true exponential dynamics. The reason 
is that, in contrast to protocol type A, the mixing pattern 
extends to the side walls of the channel (as in Fig. [ijd)) 
where the flow is less intense and stretching is slower. 
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FIG. 2. (a-c) Mixing pattern at periods 2, 3 and 6. Pat- 
terns (a) and (b) correspond to times shortly after the dye 
spot has entered the mixing region, while in (c) the dye pat- 
tern traces out the strange eigenmode. Since the flow is de- 
terministic, the supports of the escaping lobes for the three 
different times are non-intersecting domains, (d) For the same 
experiment as (a-c), blowup of a downstream lobe, where we 
have superimposed the pattern of (a) and (b) (white), and 
the long-time pattern of (c) (gray). Filaments with short 
residence times are found in white 'holes' of the long-time 
pattern. 



Therefore, dye may remain for long times in the vicin- 
ity of the walls, whereas it escapes much faster when far 
from the walls. As a result, the dye is increasingly located 
in the vicinity of the walls, and no persistent pattern is 
observed. Protocol type B is characterized by two spe- 
cific stagnation points, whose unstable manifolds divide 
the upstream flow and the mixing region (red points in 
Fig. [ijd)). These points belong to the set of periodic 
orbits, but fluid approaches them much more slowly (al- 
gebraically rather than exponentially) than other peri- 
odic points. Fluid particles are therefore exchanged more 
rapidly between the periodic orbits of the bulk than be- 
tween the wall region and the bulk, (see [39] for more 
details). This prevents the onset of a persistent strange 
eigenmode. 

In this article, we will mostly focus on the study of 
type A protocols, where the dye pattern converges to 
strange eigenmode, but we also address the case of type 
B protocols. 



IV. MIXING MEASURES DERIVED FROM 
THE EVOLUTION OF A DYE SPOT 

We now consider the different information about mix- 
ing efficiency that can be gained from our experiments, 
particularly from the strange eigenmode pattern. 



A. Location of well-mixed particles 

The eigenmode pattern can be used to give predictions 
about the efficiency of mixing for fluid particles injected 
at different locations. We illustrate this with the example 
of two dye spots represented in Fig. [sjja), for which one 
would like to know if they will be efficiently mixed by the 
rods. The degree of mixing experienced by a fluid particle 
depends on its distance to the stable manifold Wg of the 
chaotic saddle. Particles lying within a distance wb of the 
stable manifold, where wb is the Batchelor scale, are very 
well mixed since, at long times, they lie within a distance 
Wb to the unstable manifold Wu. In other words, they 
are eventually compressed onto the eigenmode pattern. 
On the other hand, particles lying further from Wg escape 
the region having experienced little mixing. 

For our stirring configurations, an approximation to 
the stable manifold Wg is easily obtained: it is the mir- 
ror image (with respect to the axis of the rods' centers) 
of the strange eigenmode pattern, which traces out the 
unstable manifold. This comes directly from (i) the time- 
reversibility of Stokes flows, and (ii) the symmetry of the 
trajectories of the rods with respect to the horizontal axis 
linking their rotation centers. Mirror images are there- 
fore to be taken when both rods lie on the symmetry axis. 
It is then possible to predict the fate of a fluid particle 
from the intersection of its mirror image (blue circles in 
Fig. [sj^a)) with the strange eigenmode. The intersecting 
part of the spot will be well mixed, as opposed to the 
non-intersecting part. In Fig. |3][a), the dye spot on the 
left is very well mixed as its mirror image is nested within 
a densely-striated lobe of the eigenmode. The dye spot 
on the right, however, escapes downstream without being 
mixed, as its mirror image has no intersection with the 
strange eigenmode. This has been checked experimen- 
tally. 

In the same way, in Fig. [sj^b), the parts of the dye spot 
whose mirror image is located on the strange eigenmode 
will be very well mixed with surrounding fluid. The other 
parts located on the mirror image of white fllaments are 
less mixed - although they may experience some stretch- 
ing - and will typically retain their initial upstream con- 
centration level once advected downstream of the rods. 
Such considerations are useful to determine the ideal lo- 
cation to inject a substance that is to be mixed. This 
ideal location, as shown above, is inside the mirror im- 
age of the strange eigenmode. 



B. Distribution of residence times 

We now consider another quantity of interest for mix- 
ing: the residence-time distribution (RTD). Such a char- 
acterization of mixing is very common in chemical engi- 
neering and was described in detail by Danckwerts |34j . 
Since the residence time of fluid particles inside the mix- 
ing region is related to the transport properties of the 
flow, and, to some extent, to its mixing properties, the 



5 






(a) 



FIG. 3. The fate of dye spots injected upstream (solid dark 
spots sketciied on (a) and (b)) can be inferred from the posi- 
tion of their mirror image across the axis of the rods (hollow 
spots), (a) Example of a type A protocol. Mirror images that 
intersect the strange eigenmode, such as the left spot, corre- 
spond to well-mixed fluid particles. The right spot, however, 
has no intersection with the strange eigenmode, and will es- 
cape downstream without being caught by the rods and with- 
out being mixed, (b) Example of a type B protocol. The 
dye spot falls on both the support of the permanent pattern 
(gray patches) and its complement (white stripes), meaning 
that only the portion located on the support of the permanent 
pattern will be smeared by diffusion. It is therefore possible to 
obtain detailed information, such as the fraction of a dye spot 
that will be well-mixed, by superimposing the mirror image 
of the spot onto the strange eigenmode pattern. 



analysis of the RTD can be used to obtain qualitative in- 
formation about mixing. We show below that an approx- 
imate quantitative measure of mixing quality can also be 
derived from the RTD by distinguishing between resi- 
dence times that are greater or smaller than the typical 
time needed for diffusion to be effective. 

In our experimental device, we can determine the RTD 
for dye particles coming from the initial dye spot by mea- 
suring the fraction of the dye that escapes the mixing 
region during the time interval from t to t + T. Gener- 
ally, the residence time distribution can be measured by 
recording the amount of dye at the outlet of the mixing 
region, either in model experiments or even in industrial 
settings. It is usually more convenient to measure the 
RTD than to obtain pictures of the strange eigenmode, 
because of obstructions or three-dimensional complexity. 

In Fig. |4|^a) we plot the RTD for a type A and B pro- 
tocols with the same rotation speed (only the rotation 
direction is different). Let us give first a qualitive analy- 
ses of the two RTD's. Two distinct regimes are observed: 
(i) a peak for short residence times, which is more pro- 
nounced for the type A protocol; followed by (ii) an ex- 
ponential decay regime (see the inset for a semilog plot 
of the RTD). For the type A protocol, most dye parti- 
cles that contribute to the short-time peak are in fact 
not entrained inside the chaotic mixing region; they flow 
almost unperturbed along the sides of the channel or, to 
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FIG. 4. (a) Experiments: Residence time distribution 
(RTD) inside the mixing region for a typical A protocol (black 
circles) and B protocol (gray triangles). A vertical dashed line 
represents idiff, the time where diffusion starts blurring fluid 
particles with different concentrations. The type A protocol 
is less effective at mixing than the B protocol, as reflected by 
the pronounced short-time peak of the RTD. Inset: RTDs on 
a semi-log scale display exponential tails, with slightly differ- 
ent slopes for each protocol type. The solid line shows the 
purely-exponential RTD of a perfectly random mixer, (b) 
Open-flow baker's map: Circles: RTD inside the mixing 
region (inset: RTD on a semi-log scale). The minimum exit 
time is three periods. Solid line: exponential RTD of a per- 
fectly random mixer with the same flow rate. 



a lesser extent, are entrained inside the mixing region 
and expelled downstream immediately. For both proto- 
cols, the onset of the exponential regime is associated 
with particles that are trapped in the mixing region for 
some time. These particles are close to the chaotic sad- 
dle, and are thus affected by their proximity to its peri- 
odic orbits [211132] (see Section III). Diffusion plays no 
role in this regime until particles have reached the res- 
idence time tdiff, the time at which dye filaments reach 
the Batchelor scale. This marks the onset of the strange 
eigenmode regime and of good mixing. 

A first quantitative index of mixing is therefore the 
fraction of fluid particles that escape only after they have 
been stretched enough to be blurred with other particles 
by diffusion: 



(2) 



with 51 the RTD. Using this index, we see in Fig. Qa) 
that the type B protocol is more efficient from a mixing 
perspective than the type A protocol: in the latter case, 
many fluid particles move along the sides the channel 
without being mixed. 
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However, this index suffers from two drawbacks. First, 
due to local disparities of the stretching rate, some dye 
filaments reach the Batchelor scale and start diffusing 
earlier than others. Specifically, some fluid particles may 
spend a long time near weakly unstable orbits and thus 
retain their initial concentration level, while most par- 
ticles are well into the diffusive regime. It is therefore 
difhcult to define a single diffusion time. Our choice for 
tdis in Fig. |4][a) is a pragmatic one: it is the period when 
most filaments start fading from pure black to intermedi- 
ate gray. The second drawback is that the RTD - hence 
Odiff - depends on the location of the injected dye spot. 
More precisely, it depends on the fraction of the spot in- 
tersecting the mirror image of the strange eigenmode, as 
was discussed above. For type A protocols, for example, 
a dye spot injected at the center of the channel yields a 
higher mixing index than a dye spot injected close to the 
channel walls. This problem can be resolved by covering 
with dye a strip that encloses all fluid particles that en- 
ter the mixing region during one period, or by dyeing all 
incoming fluid from a chosen starting time, as proposed 
in [34] . These solutions present experimental difficulties. 
In the next subsection we introduce a more refined index 
based on the strange eigenmode pattern that does not 
suffer from the two limitations mentioned above. 

In chemical engineering, distributions of residence 
times are often used to compute another mixing index, 
that is the departure from perfect mixing |34) . A perfect 
mixer is defined [U [34] as a device where an incoming 
fluid particle has a constant escape probability with time. 
It is therefore characterized by an exponential RTD, with 
escape timescale tq = Q/V, with Q the flowrate and V 
the volume of the mixing region, according to mass con- 
servation. Mixers that generate vigorous turbulence at 
high Reynolds number, for example, shuffle fluid parti- 
cles randomly and may be close to perfect mixers - al- 
though the minimum exit time is always non-zero for a 
real system. Departures from this random mixing may be 
interpreted as a bypass where fluid is whisked away from 
the mixing zone, or dead zones where fluid is trapped 
for long times. One may therefore define a mixing index 
that measures the difference between the observed and 
perfectly random residence time distributions: 



A3? : 



'Q 



-exp(-t/TQ) 



dt. 



(3) 



For type A protocols, for example, the main contribution 
to A3? comes from fluid particles advected rapidly along 
the channel sides, and for these protocols this measure 
correlates well with adiff. 

However, in the general case there is no direct relation 
between A3? and the fraction of fluid particles blurred 
by diffusion. Using A3? as a measure of mixing, it is 
even possible to find mixers that have a better mixing 
efficiency than the "perfect mixer." For example, this 
will occur if there is a minimum residence time for fluid 
particles in the mixing region - this minimum time being 
greater than the diffusion time tdiff- This is impossible in 



our experiments, since the rods transport some fluid par- 
ticles rapidly from the upstream side to the downstream 
side of the mixing region, so that very small residence 
times are unavoidable. A mixer with more rods arranged 
in an ingenious way may nevertheless increase the mini- 
mum residence time up to tdiff- 

In Fig. |4][b), we show an example of such an RTD ob- 
tained with the open-flow baker's map, where fluid parti- 
cles circulate during at least three periods before they can 
exit the mixing region (the map was introduced in |39j 
and is described in Section [VA| it is similar to an earlier 
map in |33j). As fluid particles are stretched during this 
residence time, this reduces the minimum scale of hetero- 
geneities that exit the mixer, and if tdis is small enough, 
all fluid particles will be mixed by diffusion. 

We conclude, therefore, that adiff is a more relevant 
measure of mixing quality than A3?, since the latter is not 
always correlated to the fraction of well-mixed particles. 



C. Eigenmode index 

We now present a final index of mixing quality based 
directly on statistics of the invariant concentration field 
CsE of the strange eigenmode. We show that the nor- 
malized standard deviation of Cse is a relevant measure 
of mixing quality, and we examine how this mixing index 
depends on the diffusivity of dye. 



1. Definition 

We define a new index of mixing called eigenmode in- 
dex. The eigenmode index is defined as the standard 
deviation a{CsE) of the strange eigenmode divided by 
its mean concentration {Cse)'- 



<^SE = <7{Cse)/ {Cse)- 



(4) 



asE is easily measured in decay experiments of type A 
protocols as the normalized standard deviation of the 
concentration field after a few mixing periods. This quan- 
tity becomes time-independent once the concentration 
field has achieved an eigenmode pattern, and is char- 
acteristic of the flow for a given value of diffusivity p^ . 

asE is a measure of the fluctuations of the strange 
eigenmode, in particular of the coverage of the fluid do- 
main by the support of the eigenmode. The standard 
deviation of Cse will be much greater for a few isolated 
filaments than for a more uniform coverage of space by 
the eigenmode. Small values of use therefore correspond 
to good mixing. If the strange eigenmode covers a frac- 
tion A of the fluid domain, asE can be evaluated to a 
first approximation as 



<ysE 



I- A 



A 



(5) 



which is the value asE would take for a constant value 
of Cse over the support of the eigenmode. 
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FIG. 5. (a) Eigenmode pattern observed in a type A protocol, 
(b) Same pattern blurred at a scale £ = 25 pixels, (c) Evolu- 
tion of 1 + ag^ (t) with the blurring scale £ plotted in log-log 
coordinates, where cr|g(t) is the rescaled variance of the dye 
pattern at time t. A power-law scaling is observed. According 
to Eq. ([6|, the small value of the slope (0.05) implies that the 
correlation dimension D{2) is close to the dimension of space 
do, meaning that the pattern covers the space well. The lower 
cut-off at the diffusion scale wb, as well as the upper cut-off 
imposed by the size of the largest lobes L, are shown as dotted 
lines. 



2. Dependence of the eigenmode index on the diffusivity 

An important property of the eigenmode index asE 
is its specific dependence on the diffusion coefficient, or 
equivalently, the pixel size at which the concentration 
field is probed. This dependence is directly related to 
the multifractal properties [H] of the unstable manifold 
of the chaotic saddle [111 ESI EH E^, which is the set 
that supports the strange eigenmode pattern as we saw 
in Sec. inil 

In our experiments, the temporal development of the 
mixing pattern resembles the construction of an (inho- 
mogeneous) ID Cantor set 44 . The mixing region is 
stretched and folded, and a fraction is removed and re- 
placed by dye-free fluid. This is analogous to the frac- 
tal Cantor construction where one repeatedly removes a 
fraction of each segment - the equivalent of an iteration 
of the construction step being one stirring period. The 
fractal character of the chaotic saddle and its manifolds 
has been studied previously [Ul UHl EU El] . Studies of 



fractal patterns due to chaotic advection have been per- 
formed in closed flows as well, for the classical closed-flow 
baker's map [?S1|35], the blinking vortex flow [H1I1Z]> or 
the flow between eccentric cylinders Here, we are in- 
terested in the implications of this fractal character for 
the eigenmode index. 

Let us consider the evolution oi ugE when the pattern 
is blurred at different scales. We perform the blurring 
by convolving the pattern by a constant step function of 
width I (in ID for the map and in 2D for experimen- 
tal data). This mimics the action of a diffusion coefh- 
cients k that smoothens the pattern at the diffusive scale 
£ = wb(k) = ^/k/A where A is the mean stretching (Lya- 
punov exponent). We have represented in Fig. [sja) and 
(b) the pattern of the strange eigenmode for a proto- 
col of type A, and the same pattern blurred by a square 
kernel of size 25 pixels (the size of the initial picture in 
Fig.[5ja) is 1700 x 1500 pixels). The blurring removes the 
finer details of the pattern by merging together neigh- 
boring strips. It follows from the fractal character of the 
pattern that a part of the blurred pattern of size S is 
statistically equivalent to a part of the initial pattern of 
size S/£. This property implies the following power-law 
scaling, that we derive in Appendix \K\ 



1 



2 

(^SE 



^D{2)-da 



(6) 



with do the dimension of space and D{2) the fractal cor- 
relation dimension of the support of the eigenmode. The 
quantity do^D{2) is a fractal codimension that measures 
the fraction of space that the eigenmode fails to cover. 

We check the validity of the power-law scaling ^ for 
the eigenmode pattern of Fig. [5ja) . We blur the mixing 
pattern at different scales (see Fig. [5]jb)), and compute 
the resulting intensity of segregation crgE- For the do- 
main over which asE is computed, we select only the fil- 
amentary pattern and not the large white strips close to 
the channel walls which do not participate in the stretch- 
ing and folding process. In addition, we have considered 
the concentration field inside the mixing region rather 
than in the downstream lobes, as the pattern is larger 
and has a wider range of scales so that we expect to 
see the fractal scaling on a larger range. The values of 
1 -I- cr|^ are plotted against the blurring scale £ in Fig. [5] 
for different mixing times. We observe a power-law scal- 
ing for all curves over a full 1.5 decades in space, which is 
the signature of the ongoing stretching and folding pro- 
cess. All curves converge rapidly onto a single curve as 
a result of the convergence to a persistent pattern. At 
early times the effective correlation dimension is smaller 
than the value for the persistent eigenmode, since at the 
beginning of an experiment we introduce a small spot of 
dye and the coverage of space by the dye is less important 
than for the final pattern. 

The lower limit of the fractal scaling range corresponds 
to the cut-off imposed by the physical molecular diffusion 
(as opposed to the artificial blurring we perform) in the 
experiment. In our setup, the spatial resolution of the 
camera is smaller than the diffusive scale so that blurring 



8 



FIG. 6. The action of the open-flow baker's map; g cuts the 
central square of a channel strip into two strips, stretches 
them, and re-stacks them inside the square shifted down- 
stream to mimic advection in a channel. 



the pattern below the diffusive scale does not change it. 
For large values of £, no fractal scaling is expected once 
the blurring scale is larger than the width of the largest 
hole in the pattern. We verify that this scale corresponds 
roughly to the end of the power-law scaling range. The 
scaling of 1 -f (t|^ allows one to extrapolate the value 
of the mixing index to different values of the diffusivity. 
This is of practical interest for applications where the 
substance to be mixed has a different diffusivity from 
dye used in a model experiment, or when dealing with 
effective numerical diffusivity in a numerical simulation. 



V. APPLICATION TO CONTINUOUS 
INJECTION OF DYE 



a) 
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FIG. 7. Mixing patterns obtained for the open-flow baker's 
map; (a) For a single square-shaped blob injected at f = 
and short residence times (plotted with different colors, from 
blue to red); (b) Same initial data, but for a large number 
of iterations of the map, so that the concentration pattern 
approximates the strange eigenmode. (c) Concentration for 
a continuous pulse proflle entering the mixing region. Large 
fluctuations are obtained in the same regions as in (a), that is, 
for short residence times, while fluctuations are much weaker 
on the support of the strange eigenmode visible in (b). 



We now turn to the investigation of open-flow mixing 
for the case of the continuous injection of dye. 



A. The open-flow baker's map 

Let us first introduce an idealized ID map that will 
be examined in parallel with the experiments. The map 
was described in previous studies [331 ISH] ) and serves as 
a toy model for chaotic advection in open flows. It is 
an open-flow variant of the well-studied area-preserving 
baker's map [13 HU SHI HH] ■ The action of the map is 
depicted in Fig. [6| The map is defined on a channel-like 
domain, that is an infinite strip composed of a square 
central mixing region, and corresponding upstream and 
downstream regions. At each iteration, the action of the 
map consists of two parts: (i) a translation of the strip by 
a distance U, mimicking the global advection in an open 
flow; and (ii) a 'traditional' baker's map inside the central 
mixing region. The traditional baker's map consists of a 
division into two strips, which are then compressed and 
re-stacked, preserving area. The map therefore repro- 
duces the main ingredients of chaotic advection in open 
flows: stretching, folding, and escape. 

The map has the property of mapping a distribution in- 
variant along the direction transverse to the channel into 
another such distribution. We restrict to such invariant 
distributions for the upstream concentration field, so that 
the action of the map amounts to a ID transformation g. 
Specifically, we study a family of maps defined inside the 



unit interval by 

g : X ^ gi{x)U g2{x), < x < 1, (7a) 

where 

gi{x) = U + jx; g2{x) = U + + {1 - -f)x (7b) 



and the union (U) symbol in (7a I means that g is one- 
to-two: every point x has two images given by gi{x) 
and g2{x). The parameter < 7 < 1 controls the in- 
homogeneity of stretching; 7 = 1/2 corresponding to ho- 
mogeneous stretching. Diffusion is mimicked by letting 
the concentration evolve according to the heat equation 
with diffusivity k during a unit time interval in between 
iterations of the maps. 

Such baker's maps are meant to mimic type A proto- 
cols; they cannot describe type B protocols fully because 
of the absence of stagnation points. As was observed in 
the experiments, a spot of dye is rapidly transformed by 
the map into a persistent pattern, that is, the strange 
eigenmode of the map. 



B. Case of continuous injection of dye 

For a large class of industrial processes, fluid advected 
upstream of the mixing region is always inhomogeneous, 
and its concentration field measured in a sufficiently large 
region is characterized by stationary statistics. For ex- 
ample, one could realize a mixing experiment by filling 
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one half of the upstream channel with blue dye and the 
other half with red dye, and try to mix together the blue 
and red streams inside the mixing region. 

In our experiments, we focused instead on the time- 
localized injection of a single spot of dye. Practical rea- 
sons for this choice are twofold. First, the experimen- 
tal realization of a stationary injection of dye is more 
problematic, especially if it has to be carried out until 
the downstream concentration field exhibits stationary 
statistics. Second, the choice of an upstream stationary 
field is always arbitrary: we may separate the channel 
into two blue and red halves, but we might also choose 
instead a zebra-like pattern of arbitrary wavelength. It 
is likely that any measure of mixing performed on such 
experiments will depend on some characteristics of the 
upstream pattern, for instance the typical lengthscale of 
homogeneity: a thin "premixed" zebra pattern will result 
in a more homogeneous downstream field than two paral- 
lel streams. In contrast, the strange eigenmode observed 
for a dye spot can be regarded as a unique signature of 
the underlying flow. We will now explain the relationship 
between the pattern for stationary injection of dye and 
the strange eigenmode resulting from a dye spot. 

In the case of a stationary injection of dye, the incom- 
ing fluid can be divided in distinct patches that enter the 
mixing region at different times, such as the patch con- 
taining the dye spot that we considered in Sec. |III[ Inside 
the mixing region, some of these patches have been com- 
pressed onto the support of the strange eigenmode, where 
they are blurred with other patches by diffusion. The re- 
mainder of the mixing region is covered by patches that 
have not yet been mixed with others and bear concentra- 
tion levels similar to the upstream levels. Therefore, high 
fluctuations of the concentration pattern are found on 
the complement of the strange eigenmode, whereas fluid 
is homogeneous on the strange eigenmode itself. This is 
illustrated with the open- flow map in Fig. [Tj fluctuations 
of the continuous pattern are found in the holes of the 
strange eigenmode. 



C. Intensity of segregation 

We now describe a measure of mixing commonly used 
in chemical engineering for continuous mixing, the inten- 
sity of segregation, and we examine the relation between 
this quantity and the eigenmode index. 




FIG. 8. asE vs. ac for the baker's map ([7| for different values 
of the parameters U (ranging from 4 x 10"'' to 0.2) and 7 
(ranging from 0.3 to 0.7). -I- and • symbols correspond to 
two different typical scales of the inlet pattern (0.4 and 1). A 
correlation is observed between ass and ac, becoming linear 
when they are small. The analytical relations ^ and 1 10 1 
(for poor and good mixing) are plotted as solid lines. 



(resp. downstream) region chosen to be large enough so 
that (Tc is time-independent, ac = corresponds to 
perfect mixing, and ac < 1 since mixing can only reduce 
fluctuations. This quantity was introduced by Danckw- 
erts |34j and is widely used in chemical engineering. The 
quantity 1 — ac can be regarded as the fraction of fluctua- 
tions suppressed in the mixing region, ac is independent 
of the global intensity of scalar concentration, however it 
may depend on the spatial distribution of heterogeneity 
in the upstream pattern. 



2. Link between ac and asE 

For a poor mixer with a meager coverage of space by 
the eigenmode, there are typically large holes between 
parts of the eigenmode where the greatest fluctuations of 
the continuous pattern are nested. In these holes, con- 
centration levels are of the same order as in the upstream 
pattern, since the pattern has been compressed only to a 
small extent and has not yet been smeared by diffusion. 
In such a case the value of ac is expected to be close to 
1—^1, the fraction of space occupied by the complement 
of the eigenmode. Using ([5| , we deduce the approximate 
relation between ac and asE hi the limit of poor mixing: 



1 . Definition 

For a continuous injection of inhomogeneous fluid, the 
quality of mixing can be measured by the intensity of 
segregation 



(8) 



where (Tin and (Tout are respectively the standard devia- 
tion of the concentration field measured on an upstream 



<7SE 



1 - (Tc 



for (Tc 1. 



(9) 



We now use the open-flow baker's map ([t]) to investi- 
gate the strength of the correlation between agE and ac- 
The map allows us to simulate both a continuous or lim- 
ited injection of scalar, and to also vary the parameters 
U and 7 to test different mixing hows. For the continu- 
ous injection of scalar we choose a square wave of spatial 
period A. ac is measured downstream over a domain of 
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width lOOA. We compute both ac and gse for 400 dif- 
ferent maps ([7]) where U and 7 are uniformly distributed 
in the intervals [0.02,0.2] and [0.3,0.7], respectively. The 
comparison is done for two different inlet wavelengths 
A — 0.4 and 1. The diffusivity is set to k = 10~^, which 
amounts to setting the diffusive cut-off scale I ^ 10^"*. ft 
is apparent in Fig. [8] that there is a strong correlation be- 
tween the two mixing indices, as expected from the above 
qualitative arguments. The dispersion around the mean 
curve is quite small and decreases when the wavelength 
of the inlet profile is increased. We found the root-mean- 
square deviation around the mean curve to be 8% for 
A = 0.4, and 5% for A = 1. In comparison, we obtain 
a dispersion of the order of 1% when plotting the inten- 
sity of segregation measured for the two values of the 
wavelength against each other. We find the master curve 
<^SE — f\{o'c) to depend slightly on the wavelength of 
the continuous profile A, since the intensity of segregation 
depends on A. Note that not all the results of our simu- 
lations of the map are plotted in Fig. [Sj as for the largest 
values of U, and inhomogeneous values of stretching (7 
far from 1/2), the strange eigenmode scarcely covers the 
unit interval and values of asE he outside the range rep- 
resented here. 

The evident relationship between asE and ac sug- 
gests that asE could be used as a convenient substitute 
for the intensity of segregation ac to evaluate mixing 
efficiency. Additionally, the relation between the two 
measures can be easily understood for the two extreme 
regimes of very good or extremely poor mixing. Proto- 
cols with a short mean residence time (e.g., large values 
of U for the map Q) are characterized by eigenmodes 
that cover a small fraction of space. In this limit of poor 
mixing, we expect asE arid ac to obey Eq. (|9|. The 
relation ^ is plotted in Fig. |8] for comparison with the 
data. It approximates the data well for ac ^ 0.8, that 
is, for poor mixing. 

More interesting for practical purposes is the case of 
good mixing and long mean residence times (e.g., small 
values of U for the map), for which it is important to 
accurately rate protocols according to their mixing ef- 
ficiency. The relationship between asE and ac can be 
computed for the special case of an incoming square wave, 
which has been used in the simulations. In the limit of 
large residence times, many time-periods are needed be- 
fore a full spatial half- wavelength of the inlet pattern has 
entered the mixing region. Hence, the color of injected 
fluid (i.e., black or white) is homogeneous for long times, 
as in decay experiments. The concentration pattern is 
close to the strange eigenmode for almost all times, ex- 
cept during the transition periods of the incoming fluid 
between dyed and undyed. Inside the mixing region, the 
standard deviation is therefore proportional to the mean 
concentration: 

a{C,t) ^a{C){t). 

By definition Q, the constant of proportionality a 
equals asE- The quantity {C){t) oscillates around the 



mean value of the upstream concentration on exponential 
branches corresponding to the decay of the eigenmode for 
each color. Simple algebra then leads to the relation 

asE — — for small ac, (10) 

where ctq is the standard deviation of the upstream pat- 
tern normalized by the mean concentration (ctq — 1 for 
the square waves considered here) . The above relation is 
verified on Fig. |8]for small values of ac- 



VI. CONCLUSIONS AND DISCUSSION 

In this article, we have examined quantitative criteria 
to rank open fiows with chaotic advection according to 
their mixing efficiency. Measurements performed for a 
continuous injection of dye depend strongly on the spa- 
tial distribution of the incoming dye pattern, thus ob- 
scuring the inherent efficiency of the mixing device. To 
overcome this difficulty, we have proposed a measure of 
mixing based on the strange eigenmode of the flow that 
depends only on the mixing flow and the scalar diffu- 
sivity. Our eigenmode index is defined as the standard 
deviation asE of the eigenmode and measures approxi- 
mately the area covered by fluid particles that are well 
mixed. 

As expected, we found a good correlation between asE 
and measures of the variance for the continuous injection 
of a periodic pattern with a fixed lengthscale. The mea- 
sure as E therefore ranks mixing protocols in the same or- 
der as Danckwerts' intensity of segregation. The reason 
for this correlation is that fluctuations corresponding to 
short residence times inside the mixing region are found 
on the complementary subspace of the strange eigen- 
mode. We have also noted that the dependence of asE 
on the scalar diffusivity can be predicted if the resolution 
of the eigenmode pattern is fine enough to obtain its frac- 
tal correlation dimension. As asE is easy to measure and 
has physical and practical relevance, we advocate the use 
of this new criterion for evaluating mixing efficiency. 

At first sight, asE may seem of little use for protocols 
(such as our type B protocol) where a true strange eigen- 
mode is not present because of regions of slow stretch- 
ing that trap fluid particles for long times 03]. The 
vicinity of "sticky" elliptical islands or no-slip solid walls 
are examples of regions that will be covered very slowly 
by an incoming spot of dye. Conversely, dye will be 
trapped in such regions for longer times than in the 
rest in the mixing region, so that no persistent eigen- 
mode is observed [39]. However, these regions contribute 
only weakly to fluctuations in the outflowing fluid for 
the case of a continuous injection of dye, since parti- 
cles that finally escape into the rest of the chaotic re- 
gion have typically experienced much more stretching 
than fluid particles with short residence times. Therefore, 
we suggest measuring asE at intermediate times when a 
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quasi-persistent pattern is observed, that is, when dye fil- 
aments have reached the Batchelor scale everywhere but 
on these small regions of anomalous stretching. Such a 
measure will slightly underestimate the efficiency of mix- 
ing as compared to other protocols, because some areas 
not yet covered by dye will contribute to <jse, whereas 
they do not contribute to fluctuations of the continu- 
ous pattern. Nevertheless, this is only a small effect. 
(JSE should not be measured at long times when dye is 
mostly concentrated in the slow-stretching regions, since 
the fluctuations of this pattern cannot be related to those 
for continuous injection. Most efficient mixing protocols 
display either a true strange eigenmode, or a "quasi" 
strange eigenmode corresponding to a quasi-permanent 
pattern observed during a large range of intermediate 
times. Otherwise, a large portion of the mixing region 
is occupied by fluid that is trapped for long times. This 
may happen, for example, if the rods do not come close 
enough to the channel walls for a type B protocol. As 
the mean residence time in the mixing region is deter- 
mined solely by the flowrate, this implies that there will 
be an increased fraction of short residence times, hence 
bad mixing. 

Finally, we note that our study of mixing efficiency was 
confined to a global characterization of the homogeneity 
of the outgoing fluid. In particular, no characterization of 
the spatial fluctuations of the dye pattern has yet been 
proposed for open flows. It may be relevant in prac- 
tical applications to measure the distribution of widths 
of strips of dye exiting the mixing region. Once again, 
this measure could be replaced by the simpler character- 
ization of the widths of holes in the strange eigenmode, 
where the largest fluctuations are found. This would 
require the development of dedicated image processing 
tools. 



Appendix A: Scaling of use with the difFusivity or 
box size 



Let us derive here the evolution of asE with the box 
size I, using the multifractal properties of the support of 
the strange eigenmode. 

The multifractal spectrum of a fractal set is defined 
as follows ,44 . Let p's, be the measure in the box 23 of 
size I, and g S M. The invariance of the set at differ- 
ent scales implies the following scaling of the partition 
function Z(q, £): 



(Al) 



The spectrum of multifractal dimensions is obtained from 
the spectrum T{q) as 

Dg^T{q)/iq~l). 

Using the definition of asEi we write 



1 + acp^ — „ 

(E.G.)' 



(A2) 



with Ci = C(xi) the concentration level in pixel i. To 
find the multifractal scaling, instead of summing over 
pixels we sum the concentration over larger square boxes 
23 of size i, over which the concentration has been blurred 
and is nearly a constant C'('B). Since 

with do the dimension of space and g G M, we obtain 
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with the measure p's, = C/^^ C. Using the definition 
of the multifractal spectrum, we finally get 
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